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Highlights: 

e Three Alpha-CoV strains were fully sequenced by NGS method 

e The Italian strains were classified into two novel Alpha-CoV species 

e The phylogenetic analysis on RdRp fragment sequences showed correlation to European 


strains 


Summary 


Coronaviruses (CoVs) have been detected worldwide in several bat species, which are considered the main reservoir. 
The attention to the high diversity of CoVs hosted by bats has increased during the last decade due to the high number 
of human infections caused by two zoonotic Beta-CoVs, SARS-CoV and MERS-CoV, that cause several respiratory 
diseases. Among coronaviruses, two Alpha-CoV strains (HuCoV-229E and HuCoV-NL63) cause mild respiratory 


disease that can change to severe disease in children, elderly and individuals affected by illnesses. Phylogenetic analysis 
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conducted on bat Alpha-CoV strains revealed their evolutive correlation to human strains, suggesting their origin in 
bats. The genome of CoVs is characterized by a high frequency of mutations and recombination events, increasing their 
ability to switch hosts and their zoonotic potential. In this study, three strains of Alpha-CoV genera detected in Italian 
bats (Pipistrellus kuhlii) were fully sequenced by Next Generation Sequencing (NGS) and characterized. The complete 
genome analysis showed the correlation of the Italians strains with a Chinese strain detected in 2013 and, based on CoV 
molecular species demarcation, two new Alpha-CoV species were established. The analysis of a fragment of the RNA- 
dependent RNA polymerase (RdRp) showed the correlation of the Italian strains with CoVs was only detected in the bat 


Pipistrellus genera (Pipistrellus kuhlii and Pipistrellus Pipistrellus) in European countries. 
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1. Introduction 


Bats are considered the natural reservoirs of several emerging and re-emerging viruses, such as Nipah virus, Marburg 
virus, rabies virus and coronaviruses, that have caused outbreaks in both humans and animals (Shi, 2013; Smith and 
Wang, 2013). The ecological features of bats, including their ability to fly long distances, their longevity, their large 
social colonies and their potential interactions with humans or livestock animals, facilitate virus maintenance and 
transmission, increasing the risk of intraspecies or interspecies jumping (Calisher et al., 2006). Among bat viruses, in 
the last decade, a large diversity of coronaviruses has been detected, exceeding the diversity seen in other mammalian 
hosts (Drexler et al., 2014). Coronaviruses (CoVs) (order Nidovirales, family Coronaviridae, subfamily Coronavirinae) 
are enveloped viruses characterized by a positive-sense single-stranded RNA genome of approximately 26 to 32 
kilobases and classified into four genera (Weiss and Leibowitz, 2011). Alphacoronavirus (Alpha-CoV) and 
Betacoronavirus (Beta-CoV) infect several mammal species, including humans, bats and pigs, while 
Gammacoronavirus (Gamma-CoV) and Deltacoronavirus (Delta-CoV) infect birds, wild felines, pigs and some marine 
mammal species (Woo et al., 2009b; Woo et al., 2012). The genome of coronaviruses is characterized by high 
frequency recombination and a high mutation rate, which increases their potential for interspecies and intraspecies 
jumping (Lai, 1992; Holmes, 2009). Six CoV strains are recognized to infect humans. Two Alphacoronaviruses 
(HuCov-229E, -NL63) and two Betacoronaviruses (HuCoV-OC43, -HKU1) are responsible for the common cold and 
severe respiratory pathologies in infants, elderly people and immunocompromised patients and are characterized by 
human-to-human transmission (Hu et al., 2015). The other two Betacoronaviruses species, the Severe Acute Respiratory 
Syndrome virus (SARS-CoV in 2002-2003) and the Middle East Respiratory Syndrome virus (MERS-CoV in 2012) 
caused severe respiratory pathologies with case fatality rates of 9% and 35%, respectively (WHO, www.who. int). 


Phylogenetic analysis on strains detected in bats, humans and other mammals suggested that the origin of these CoVs 
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was in bats. The Rhinolophus bat species are considered the main reservoir for SARS-related CoVs. Bat MERS-related 
CoVs were also detected in African, Chinese and Italian bats (Annan et al., 2013; Ithete et al., 2013; Lau et al., 2013; 
Corman et al., 2014; Moreno et al., 2017), supporting the hypothesis of the bat origin. In three recent studies, related 
strains of the HuCoV-229E were detected in Hipposideros bats and strains of HuCoV-NL63 were detected in the 
American tricoloured bat (Perimyotis subflavus) and Kenyan Triaenops afer species, suggesting bats as potential 
reservoir host of Alphacoronavirus human strains (Pfefferle et al., 2009; Huynh et al., 2012; Corman et al., 2015). In 
addition, relatives of HuCoV-NL63 can be grown in immortalized bat cell lines, suggesting their potential association 
with bats (Huynh et al., 2012). This has led to speculations about an evolutionary origin of all mammalian CoVs in bat 
hosts (Woo et al., 2009a; Woo et al., 2009c). However, how humans become exposed to remote wildlife viruses is not 
always clear (Wolfe et al., 2007). 

In Europe, several studies described the presence of CoVs in bat populations detecting both Alpha-CoVs and Beta- 
CoVs in Germany, Spain, Luxembourg, Italy, The Netherlands, the United Kingdom, France and Hungary (Gloza- 
Rausch et al., 2008; Reusken et al., 2010; Falcon et al., 2011; August et al., 2012; Lelli et al., 2013; Kemenesi et al., 
2014; Goffard et al., 2015; Monchatre-Leroy et al., 2017; Pauly et al., 2017) from more than 20 different bat species. 
The detection of the same CoV strains (100% nucleotide identity) in different colonies of the same bat species or the 
circulation of different genera of CoVs (Alpha-CoVs and Beta-CoVs) in the same bat species confirm the high 
heterogeneity of CoVs in bats and that bat-CoV diversity depends more on the species-specificity than the geography 
and sampling location. 

However, these studies were based on the analysis of a fragment of the RNA-dependent RNA polymerase (RdRp) gene 
that allows the assignment of the strains to the genera and not to the species. The International Committee on 
Taxonomy of Viruses (ICTV) established a molecular demarcation method for species assignment using the conserved 
domains of replicase polyprotein and the pairwise amino acid distance of 9096 as threshold value. The Alpha-CoVs are 
classified into 11 species, 6 of which detected in bats: Miniopterus bat coronavirus 1, Bat coronavirus CDPHEI5, 
Miniopterus bat coronavirus HKUS, Rhinolophus bat coronavirus HKU2, Bat coronavirus HKUIO, and Scotophilus bat 
coronavirus 512, and some strains that to date are not assigned. However, the number of bat species that host CoVs is 
still unknown and increases proportionally with the increasing of surveillance. In this study, we describe the full 
genome sequencing by Next Generation Sequencing (NGS), the characterization and the classification of two novel 


Alpha-CoV species detected from three Italian Pipistrellus kuhlii bats (Lelli et al., 2013). 


2. Materials and methods 


2.1 Sampling 


Two bat faecal samples and one carcass from three bat Pipistrellus kuhlii species were provided by a rehabilitation 
centre from Northern Italy between 2010 and 2015 and the bats species were identified according to the European bat 
identification keys based on their morphologic characteristics (Dietz, 2013). Faecal and organ samples positive for 


Alpha-CoV genera by a pan-coronavirus one-step RT-PCR (Lelli et al., 2013) were chosen for NGS analysis. 


2.2 Whole-genome sequencing 

Libraries were prepared following the sequence independent single primer amplification method (SISPA) (Djikeng et 
al., 2008). The RNA, extracted as previously described by Lelli et al. (2013), was retro-transcribed using the 
SuperScript IV Reverse Transcriptase (Invitrogen, Monza, Italy), starting with 9 ul of RNA and following the 
manufacturer's instructions. Twenty microlitres of cDNA were used to synthesize the second strand of cDNA by DNA 
Polymerase I Large (Klenow) Fragment (Promega, Milan, Italy) and then amplified by the Expand High Fidelity PCR 
System (Sigma Aldrich S.R.L., Milan, Italy). The PCR amplicons were purified using one volume of Agencourt 
AMPure XP beads (Beckman, Milan, Italy) following the manufacturer's instructions and eluted in 40 ul of nuclease- 
free water. Five hundred nanograms of purified DNA, quantified with the Qubit dsDNA HS Assay Kit (Thermo Fisher 
Scientific, Rodano, MI, Italy) were digested with the EcoRV enzyme (New England BioLabs, Pero, MI, Italy) and then 
purified with a 1.8x volume of Agencourt AMPure XP beads (Beckman, Milan, Italy). The libraries were prepared by 
NEBNext Fast DNA Library Prep Set for Ion Torrent following the standard protocol for 100 ng of DNA. The barcoded 
libraries were mixed, and the pool was used for the Emulsion PCR performed by the Ion PGM Hi-Q OT2 Kit. The 
sequencing run was performed according to the manufacturer's instructions (Ion PGM Hi-Q Sequencing Kit) (Thermo 


Fisher Scientific) by Ion Personal Genome Machine (PGM) on the Ion 318 Chip v2. 


2.3 Genome structure and phylogenetic analyses. 

NGS, previously described by Moreno et al. (2017) was applied in order to obtain the complete genome. Data obtained 
by the Ion Torrent sequencer were analysed by the online portal Galaxy Aries (https://aries.iss.it). The reads were 
checked, cleaned up, and trimmed, and the sequences shorter than 50 nt were filtered. Host sequences were removed by 
mapping the reads against the Megabat and Microbat complete genomes downloaded from Genome Browser 
(https://www.genome.ucsc.edu) using the Bowtie 2 tool. The reads aligned by the BLASTn tool to the bacterial non- 
redundant nucleotide database RefSeq (https://www.ncbi.nlm.nih.gov/refseq/; E-value »10 5) were removed by Galaxy 
Aries. The remaining reads were aligned with the viral non-redundant nucleotide database RefSeq 
(https://www.ncbi.nlm.nih.gov/refseq) and were parsed with the MEGAN6 software. The reads that showed no 


significant hits to the reference database were assigned to the unclassified reads. 


The sequences classified into the Coronaviridae family were extracted and assembled into contigs by a de novo 
assembling method, using the default parameters, and excluding those shorter than 1000 bases using SPAdes tool 
(Galaxy Aries). The closest viral sequences were chosen by a BLASTn analysis and used to map the reads by the online 
tool Bowtie2 (Galaxy Aries). The output was visualized by the Integrative Genomics Viewer (IGV) software 
(http://software.broadinstitute.org/software/igv/), and the consensus sequence was extracted. Nucleotide and amino 
acids sequences were aligned, and the pairwise identity values were calculated with MEGA7 software 
(www.megasoftware.net). The open reading frames (ORFs) were predicted using the online tool ORF Finder (NCBI, 
http://www.ncbi.nlm.nih.gov/gorf/gorf.html). The potential cleavage sites in the orflab polyprotein were predicted by 
amino acid sequence alignment with other CoV strains and by using the online tool NetCorona 1.0 Server 
(http://www.cbs.dtu.dk/services/NetCorona/) (Kiemer et al., 2004). Comparison of the sequence distances of BatCoV- 
Ita4 and the closest Alpha-CoV sequences were confirmed using SSE v1.2 (Simmonds, 2012). 

A dataset of complete genome references for Alpha-CoVs and Beta-CoVs species from the ICTV taxonomy report 
(https://talk.ictvonline.org) were obtained including the full genome sequences that displayed the highest nucleotide 
similarity to strains sequenced in this study, which resulted in 59 CoV sequences. A second dataset was used to build a 
ML tree using the partial sequence of the RdRp gene (409 nt) sequenced worldwide, excluding identical strains from 
the same study and bat species and resulting in 226 CoVs sequences. In both ML trees Beta-CoVs from different 
species were used as an outgroup. 

To test the presence of recombination by RDP4 (Martin et al., 2015), six different methods were applied: GENECONV, 
BootScan, MaxChi, Chimaera, 3Seq, and SiScan, using the default settings. 

The Maximum likelihood (ML) phylogenetic trees were built using MEGA7 software, applying, as a substitution 
model, a general time-reversible (GTR) model with a gamma-distributed (G) rate variation across sites, a proportion of 
invariant sites (I) (GTR+G+D) and a bootstrap analyses of 1000 pseudo-replicates. 

The Bayesian phylogenetic trees were carried out using MrBayes ver. 3.1.2 (Huelsenbeck and Ronquist, 2001) using the 
sequences of the predicted proteins and excluding the most divergent strains. The Metropolis-coupled Markov chain 
Monte Carlo (MCMC) was used, starting from a random tree, run for 500 thousand heuristic search generations, 
sampling every 1,000 generations and discarding 25% of the samples as burn-in. 

Analysis of the protein families of spike proteins and the prediction of the secondary structure were performed by the 
online tools: PFAM, InterProScan, the TMHMM program (http://www.cbs.dtu.dk/services/sTMHMM/) (Apweiler et al., 
2001; Bateman et al., 2002), Phyre2 (http://www.sbg.bio.ic.ac.uk/phyre2/html/page.cgi?id=index) and the Swiss model 


(https://swissmodel.expasy.org). 


3. Results 

3.1 NGS data analysis 

The NGS run produced approximately 4 million reads. A total of 1,376,444 reads were obtained for Bat- 

CoV/P. kuhlii/Italy/206645-4 1/2011 (BatCoV-Ita3). Of these reads, 8% were host reads, 78% were bacterial reads and 
84,069 (6%) were viral reads, of which 81,069 (5.8%) were classified as Coronaviridae and 8% were unclassified. For 
Bat-CoV/P. kuhlii/Italy/3398-19/2015 (BatCoV-Ita4), 970,190 reads were retrieved, of which 6.8% were host reads, 
38% were bacterial, 72,127 (7.4%) were viral, of which 5,484 (0.5%) were Coronaviridae, and 52% were unclassified. 
For Bat-CoV/P.kuhlii/Italy/206679-3/2010 (BatCoV-Ita5), 1,602,274 reads were obtained. Of these, 0.3% were host 
reads, 98% from bacteria, and 22,621 were viral (1,496), of which 19,053 (1.1%) were Coronaviridae sequences and 
0.3% were unclassified. 

The reads classified into the Coronaviridae family were used to assemble the contigs obtaining | contig of 
approximately 27,000 nt for BatCoV-Ita3, 6 contigs > 4,000 nt for BatCoV-Ita4 and BatCoV-Ita5 contigs 22,800 nt for 
BatCoV-Ita5. The three assembled full genomes showed an average coverage of 751x, 50x and 181x, for BatCoV-Ita3, 
BatCoV-Ita4, and BatCoV-Ita5, respectively. The RdRp sequence of BatCoV-Ita4 was not obtained by the Sanger 
method used by Lelli et al. (2013). The other two Italian bat RdRp sequences showed 99% nucleotide identity with the 


RdRp region of the complete genomes obtained by NGS. 


3.2 Genome organization 


The complete genome sizes were 27,862 nt for BatCoV-Ita3, 28,129 nt for BatCoV-Ita4 and 28,146 BatCoV-Ita5, with 
a G+C content of 42%, 40.3% and 40.4%, respectively. The first sequence analysis was performed by BLASTn, 
comparing the Italian strains with those available online (https://www.ncbi.nlm.nih.gov). The BLASTn search showed 
similarities with an unclassified strain BENv-AlphaCoV/SC2013 (KJ473809), and with those viruses classified into 
HKU1O0 bat and Porcine epidemic diarrhoea virus (PEDV) species. BatCoV-Ita4 and BatCoV-Ita5 shared a 97% 
nucleotide identity (nt. id.), and BatCoV-Ita3 shared 70.9% and 71% nt. id. with BatCoV-Ita4 and BatCoV-Ita5, 
respectively, at the full genome level. The same differences were observed when the ORF nucleotide sequences were 
aligned separately. BatCoV-Ita4 and BatCoV-Ita5 shared >97% nt. id. in all the ORFs. BatCoV-Ita3 showed the highest 
differences in the S, ORF3 and N genes with <65% nt. id.. Fewer differences were observed at the ORF1ab, M and E 
genes (27046 nt. id.) compared to the other two Italian strains. Their genome organization was similar to other Alpha- 
CoV species, comprehending 6 ORFs and two non-translated termini in the order of 5’ terminus-ORFl1ab-spike-ORF3- 
envelope (E)-membrane (M)-nucleocapsid (N)-3' terminus (Table 1). In the ORFlab, it has been observed that the 


predicted slippery sequence *UUUAAAC'" is involved in the synthesis of the replicase pplab polyprotein by ribosomal 


frameshift, a characteristic of the Nidovirales order. The sizes, the genomic localization and the 15 expected cleavage 
sites of the nonstructural protein (NSP 1-16) that are encoded by ORF lab, were predicted by sequence comparison with 
other Alpha-CoV species (Table 2). BatCoV-Ita4 and BatCoV-Ita5 showed the same sequences of cleavage sites. 
BatCoV-Ita3, compared to the other 2 strains, showed two amino acid changes in the cleavage sites between 
NSPI/NSP2 and NSP12/NSP13. A leader predicted transcription regulatory sequences (TRS-L), and the putative body 
TRSs, representing signals for the discontinuous transcription of subgenomic mRNAs (sgmRNAs), have been identified 
in the three genomes (Table 1). The TRS-L and TRSs preceded the codon start of all ORFs in BatCoV-Ita3 and 
suggested the synthesis of 6 monocistronic subgenomic mRNAs. The lack of TRSs before the ORF3 gene codon start in 
BatCoV-Ita4 and BatCoV-Ita5 suggests the synthesis of 4 monocistronic and 1 polycistronic subgenomic mRNAs. The 
differences at the nucleotide level were also confirmed at the amino acid level. BatCoV-Ita4 and BatCoV-Ita5 showed 
high similarities (<97%) and high differences with BatCoV-Ita3 in the spike, ORF3 and nucleocapsid proteins. The 
ICTV has established the 90% amino acid sequence identity of the seven concatenated domains within the ORFlab as 
the threshold value to assign two strains to the same species: NSP3 (ADRP), NSP5 (3CLpro), NSP12 (RdRp), NSP13 
(Hel, NTPase), NSP14 (ExoN, NMT), NSP15 (NendoU), and NSP16 (OMT). To classify the Italian strains into known 
coronavirus species, the ORFlab concatenated domains were compared with the 11 Alpha-CoV species: Miniopterus 
bat coronavirus 1, Bat coronavirus CDPHE15, Miniopterus bat coronavirus HKUS, Rhinolophus bat coronavirus 
HKU2, Bat coronavirus HKU10, Scotophilus bat coronavirus 512, PEDV, HuCoV-229E, HuCoV-NL63, and 
AlphaCoVs1, and some strains that to date are not assigned. BatCoV-Ita3 concatenated domains showed sequence 
identities «83.846 with all the Alpha-CoV strains. BatCoV-Ita4 and BatCoV-Ita5 shared 99.3% identity and had <79.1% 
with all other Alpha-CoVs, suggesting that the classification of the Italian strains should be into two novel Alpha-CoVs 


species. 


3.3 Phylogenetic analyses 

The RDP4 recombination detection methods, applied to the dataset of CoVs complete genomes to detect the occurrence 
of recombination, supported the absence of recombination between the Italian strains and the Alpha-CoVs strains (P 
values >0.05). As shown in the ML tree built with complete genomes (Fig. 1), the Italian strains clustered with the 
Chinese strain BtNv-AlphaCoV/SC2013 (KJ473809) out of the monophyletic clade formed by the complete genomes of 
the HKU-8, 1A, 1B and HKU10 species. The former cluster is divided into two sub-clusters: one sub-cluster 
represented by BatCoV-Ita3 and BtNv-AlphaCoV/SC2013, sharing 75% nt. id., and the other sub-cluster represented by 
BatCoV-Ita4 and BatCoV-Ita5, sharing 71% nt. id. with the Chinese strain. 

The Italian strains showed approximately 62% nt. id. with the strains classified into the HKU10 species (Hipposideros 


bat coronavirus HKU10 isolate LSH5A, Rousettus bat coronavirus HKU10 isolate 183A) and «6046 with all other 


AlphaCoV strains (Supplementary Fig. 1, Supplementary Fig. 2). Additionally, at the amino acid level, the Italian 
strains showed the highest identities with the Chinese BtNv-AlphaCoV/SC2013 strain with respect to the other Alpha- 
CoVs. BatCoV-Ita3 showed high identities in all predicted proteins excepting in the ORF3 and N proteins. BatCoV-Ita4 
and BatCoV-Ita5 showed lower identities with respect to BatCoV-Ita3, which showed high identities in the orflab 
polyprotein and M proteins (>75%) and low identities in the other predicted proteins. 

The Bayesian trees, built using the predicted protein sequences of E, M and N, confirmed the clustering of the Italian 
strains with the Chinese strain BtNv-AlphaCoV/SC2013 (data not shown). The tree built with S protein sequences 
showed a uniquely supported clade, containing the Italian strains, the BtNv-AlphaCoV/SC2013, HKU10, 1A, 1B, and 
HKUS species and the unclassified strain BtMr-AlphaCoV/SAX2011, suggesting correlation only between those bat 
species (Supplementary Fig. 3). 

The Italian bat strains showed low identities with the HuCoV-229E (<49%) and HuCoV-NL63 (<45%) strains at the 
spike protein level and had <45% identity with HuCoV-229E and <35% with HuCoV-NL63 at the Receptor Binding 
Domain (RBD) level. The prediction structure of the spike protein showed a type I membrane glycoprotein divided into 
two subunits (S1 and S2), as other Alpha-CoVs Spike proteins with most of the protein exposed on the outside of the 
virus and two transmembrane domains located at the C terminus. However, the Italian strains did not exhibit significant 
or supported similarities to the known secondary structure receptor-binding domains (HuCoV-229E, -NL63) using the 
online tool Phyre2 or the Swiss model due to their high divergences (data not shown). 

To investigate the correlation among Alpha-CoV strains previously detected worldwide, a phylogenetic tree of the 
partial RdRp gene was built (Supplementary Fig. 4). The ML showed that strains detected in the same continent shared 
>89% nt. id. and were correlated, forming monophyletic clusters while sequences from a different cluster showed a nt. 
id. «8596. Most of the Alpha-CoVs species were detected in the same continent as the 1A, 1B, Bat-CoV 512, HKU2, 
and HKUS species in Asia or the CDPHE15 species in North America. Bat coronaviruses related to human HuCoV- 
229E were retrieved in Africa and the coronaviruses related to HuCoV-NL63, in Africa and America. The HKU10 CoV 
strains showed sequences similar to those detected in Asia and Europe. Some strains formed a cluster outside of those 
classified into known Alpha-CoV species. 

The Italian strains formed two clusters with the Chinese strain BtNv-AlphaCoV/SC2013 and some European strains. At 
the RdRp partial gene level, the BtNv-AlphaCoV/SC2013 strain showed 83% nucleotide identity with BatCoV -Ita3, 
83.8% with BatCoV-Ita4 and 82.8% with BatCoV-Ita5. The first cluster is formed by BatCoV-Ita4 and BatCoV-Ita5, 
one Italian strain and one Spanish strain (P.kuh/Iprima/Spain/2007, HQ184058), collected from the bat Pipistrellus 
kuhlii species in the Southwest Piedmont region in Northern Italy (Pkuh605, KY780383) and in Spain in 2014 and 


2007. These strains shared >96.7% nt. id. The second cluster contains the BatCoV-Ita3 with one Italian strain collected 


in the centre of the Piedmont region in Northern Italy (Ppip1015C, KY780385), and a French strain (KT345294, 

Pipl Cr FR 2014), both collected in 2014. Those strains formed a monophyletic clade with two European strains, 
detected in Bulgaria (GU190239, BNM98-30/BGR/2008) and Spain (HQ184057, M.myo/I/Spain/2007), and two strains 
from South Africa (KF843855, BtCoV/GrNC1/Neo; KF843862, BtCoV/GrNC8/Neo) from the Nyctalus leisleri, Myotis 


myotis and Neoromicia capensis species, sharing with BatCoV-Ita3 approximately 83% nt. id. 


4. Discussion 


In this study, three Alpha-CoV strains from the Pipistrellus kuhlii bat species were fully sequenced. The P. kuhlii 
species is one of the most frequently described bat species in Italy that forages in urban and agricultural areas (Russo 
and Jones, 2003; Ancillotto et al., 2016). 

To fully characterize the three Alpha-CoV strains, the NGS method previously described by Moreno et al. (2017) was 
applied successfully, obtaining the Alpha-CoV complete genome sequences with high coverage rates. However, the 
lack of European Alpha-CoV complete genomes make difficult to conduct a comprehensive genetic and phylogenetic 
analysis. The analysis on the full BatCoV-Ita sequences showed similarities to the Alpha-coronavirus genera and 
genome organization with 6 open reading frames (ORFs) and the 5’ and 3’ non-translated sequences. The phylogenetic 
analysis using the complete genomes showed correlation but with a low nucleotide identity with a Chinese strain 
detected in 2013 in the Nyctalus velutinus species. The phylogenetic analysis on amino acidic sequences also confirmed 
the correlation with the Chinese strain and supported the hypothesis that bat strains of Miniopterus bat coronavirus 1, 
Miniopterus bat coronavirus HKUS, and Bat coronavirus HKU10 species and some unclassified strains may share a 
common spike ancestor. However, the analysis of the protein structure was hampered by the lack of similar spike 
protein structure. Indeed, due to the high genetic divergences with human strains it was impossible to predict the spike 
structure and the affinity with the human receptor. 

The ICTV has established that viruses sharing more than 90% amino acid sequence identity in the conserved 
concatenated domains of the orflab polyprotein can be assigned to the same CoV species 
(https://talk.ictvonline.org/taxonomy/). The ICTV demarcation criteria for genera and species allowed us to classify the 
BatCoV-Ita into two novel Alpha-CoVs species. Our results support previous findings about the high heterogeneity of 
CoVs hosted by bats and support the idea that novel species may be found in the future with increasing surveillance. 
Several studies described the presence of Alpha-CoV and Beta-CoVs in bats worldwide (Falcon et al., 2011; Gouilh et 
al., 2011; August et al., 2012; Goffard et al., 2015; Asano et al., 2016; Fischer et al., 2016; Goes et al., 2016; Subudhi et 
al., 2017; Ar Gouilh et al., 2018; Geldenhuys et al., 2018). However, most of these studies reported phylogenetic 
analysis on short sequences within the RdRp region, establishing the correlation with other CoV strains but not the 


assignment to CoV species as established by ICTV. 


In Europe, CoV strains were detected in samples from more than 20 different bat species. In Italy, a large variety of 
CoV strains were detected in Myotis nattereri, Myotis daubentonii, Myotis myotis, Rhinolophus hipposideros, Hypsugo 
savii, Pipistrellus kuhlii, Pipistrellus pipistrellus, Nyctalus noctula, Epseticus serotinus, Myotis blythii, Myotis 
oxygnathus, and Plecotus auritus (Lelli et al., 2013; De Benedictis et al., 2014; Rizzo et al., 2017). 

The phylogenetic analysis on the RdRp region showed the correlation of the Alpha-CoV strains detected in the same 
continent. Interestingly, within each geographic area most of the strains hosted by the same bat genera cluster together, 
confirming the CoVs -host coevolution. BatCoV-Ita4 and BatCoV-Ita5 strains showed high nucleotide identity with one 
Italian strain and one Spanish strain detected in 2014 and 2007, respectively, from the bat Pipistrellus kuhlii species 
(HQ184058.1; KY780383.1) (Falcon et al., 2011; Rizzo et al., 2017). The BatCoV-Ita3 result correlated with one 
Italian and one French strain (KY780385.1; KT345294.1), both collected in 2014, from the bat Pipistrellus Pipistrellus 
species (Goffard et al., 2015; Rizzo et al., 2017). The high identities at the RdRp gene level and the clustering of the 
European strains with the Italian strains suggest that the two novel Alpha-CoV species detected in this study may infect 
at least two bat species of the Pipistrellus genera (Pipistrellus kuhlii and Pipistrellus Pipistrellus) from different 
European countries. 

In contrast, some geographical clusters were represented by strains detected in different bat genera, attesting to the 
capability of the CoVs interspecies jumping that may occur when different species of bats share same roost (Leopardi et 
al., 2018). 

Bat behaviour, including flying long distances, living in large colonies, having social interactions, and cohabitating with 
different bat species, favour the interspecies or intraspecies transmission of viruses(Calisher et al., 2006). During the 
last fifteen years, two Beta-CoVs, SARS-CoV and MERS-CoV, have jumped from bats to a mammalian intermediate 
host to humans (Field, 2009; Omrani et al., 2015). In addition, strains related to human Alpha-CoVs (HuCoV-229E, 
HuCoV-NL63) have been detected in bats, indicating the importance of the bat as a CoV reservoir (Pfefferle et al., 
2009; Corman et al., 2015; Tao et al., 2017; Waruhiu et al., 2017). In this study, we characterized two Alpha-CoVs 
from Italian bats divergent from human CoVs strains and two new Alpha-CoV species. In addition, the RdRp 
phylogenetic tree showed that the strains here described were not related to the Alpha-CoV species established so far. 
This result highlights that the heterogeneity of CoVs in the bat may be higher than what is known to date. Indeed, to 
better understand the CoV species circulating in bats, their evolution and our understanding of the mechanisms 
important to cross the species barrier, it is important to have long-term vigilance followed by the complete genome 


characterization. 
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Fig. 1. Maximum phylogenetic tree based on 47 Alpha-CoVs and 12 Beta-CoVs complete genomes. The tree 


was inferred under the GTR + G + I substitution model and 1000 bootstrap resampling 


process replications showing values > 70. BatCoV-Ita3, BatCoV-Ita4 and BatCoV-Ita5 are 
reported in bold and the sequences can be retrieved under accession numbers (MH938448, 


MH938449, MH938450) 


KJ473799 BtMEAlphaCoV/F 12012 
vc 4473800 BtMEAlphaCoV/HeN2013 
KJ473798 BtMF-AlphaCoV/HuB2013 HKUS uf 
KJ473797 BIMFAlphaCoV/GD2012 
10L NC 010438 Bat coronavirus HKUB 
EU420137 Bat coronavirus 18 
Ku 3796 BtMEAlphacovinzoi2 NA T 
sp KJ473795 BtMFAlphacoviAH2031  |1B 
NC. 010437 Bat coronawrus 1A 
soor 2989269. 1 Hipposideros bat coronavirus HKU10 isolate LSHSA 
JQ989270.1 Rousettus bat coronavirus HKU10 isolate 183A | HKU10 a 
oof KJ473807.1 BIRFAlphaCoV/HuB2013 
KJ473810,1 BtMs-AlphaCoV/GS2013 
K473809.1 BINvAlphaCoV/SC2013 
BatCoV-Ita3 
BatCoV-Ita4 
BatCoV-ItaS 
AF353511 Porcine epidemic diarrhea virus strain PEDV ] PEDV - 
NC. 009657 Scotophilus bat coronawrus 512 1512 af 


































100 [ KY799179 Myotis lucifugus coronavirus Jeophe 


NC 022103 Bat coronavirus CDPHE 15/USA/2006 
KJ473806 BtMr-AlphaCoV/SAX2011 
AY567487 Human Coronavirus NL63 
KY073744 NL63-elated bat coronavirus strain BIKYNL63-9a Jus za” Ñ 
KY073745 NL63-related bat coronavirus strain BIKYNL63-9b 
KY073747 229E-related bat coronavirus strain BIKY229E-1 
KT253269 229E-related bat coronavirus isolate BtCoV/KW2E-F151/Hip cf. rub/GHA/2011 " P 
NC. 002645 Human coronavirus 229E 
KU291449 Camel aiphacoronavirus Camel229E isolate Camel229E-CoV/KCSP/KENI2015 | o 
JO410000 Alpaca respiratory coronavirus isolate CAO8-1/2008 
KT253324 Camel alphacoronavirus Camel229E isolate Camel229E -CoV/JC50/KSA/2014 kud 
MF593473 Camel alphacoronavirus Abu Dhabi B38 
NC. 028752 Camel alphacoronavirus isolate camel/Riyadh/Ry 141/2015 
oof 473808 BIREAlphaCoV/YN2012 
MF370205 Rhinolophus bat coronavirus HKU2 | HKU2 Af 
NC 009988 Bat coronavirus HKU2 
KM347965 Ferret coronavirus isolate FRCoV 
NC_023760 Mink coronavirus strain MCoV 
AY994055 Feline infectious peritonitis virus FIP 


(00848676 Feline coronavirus strain FCoV Alphacoronavirus 1 
KP849472.1 Alphacoronavrus 1 





DQ811787 Porcine Respiratory Corona Virus Infection PRCV 
NC_035191 Wencheng Sm shrew coronavirus isolate Xingguo-101 

NC_032730 Lucheng Rn rat coronavirus isolate Lucheng-19 

NC. 034972 Coronavirus AcCoV-JC34 

EF065505.1 Bat coronavirus HKU4-1 

EF065509.1 Bat coronavirus HKUS-1 

JXB69059 Middle East Respiratory Syndrome MERS 

AY278741 Severe acute respiratory syndrome SARS 

EFO65513.1 Bat coronavirus HKUS-1 

AY597011.2 Human coronavirus HKU1 Beta-CoVs 

AC 000192.1 Murine hepatitis virus strain JHM 

EF446615.1 Equine coronavirus ECoV 

DQ011855.1 Porcine hemagglutinating encephalomyelitis virus PHEV 

voc|p AY585228.1 Human coronavirus OC43 

MF593476, 1 Camel coronavirus HKU23 

U00735.2 Bovine coronavirus BCoV 
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Table 1 Locations of predicted ORFs, protein seguences, putative leader TRS-L and TRS-B 





ORF Caa No. of amino acids Sequence? 

BatCoV-Ita3 

ORF lab (TRS-L) 296-20166 6623 00067C LA A A Coooz3 

Spike 20153-24274 1373 20100— —— T —20106 

ORF3 24274-24924 216 24230— — —— 24236 

E 24947-25174 75 24939— — — 24945 
25180-25863 227 25170— — — 25176 

N 25870-27174 434 25859— —— T 25865 

BatCoV-Ita4 

ORF lab (TRS-L) 281-20175 6627 ooost CTAAAC 0057 

Spike 20172-24371 1399 20161--C———20167 

ORF3 24371-25024 217 

E 25103-25330 75 25077— —— T —25083 
25343-26050 235 25306 —— — — 25332 

N 26058-27356 432 26046— — — — 26052 

BatCoV-Ita5 

ORF 1ab (TRS-L) 253-20192 6645 0023 CTAAACoo029 

Spike 20189-24388 1399 201783--C———20184 

ORF3 24388-25041 217 

E 25120-25347 75 25094— —— T —25100 

M 25360-26067 235 25343—————5 25349 

N 26075-27313 432 26063— — —— 26069 





è Dashes represent identical nucleotides compared to the leader TRS 
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Table 2 Prediction of the putative pp lab cleavage sites 








NSP BatCoV-ITA3  BatCoV-ITA4  BatCoV-ITA5 Putative functional domain(s)* 
NSPI M!-A M'-G!? MI-G? 

NSP2  GoÓ.gG71 G8. 571 G5. G71 

NSP3 G2. 52961 G7?.G2374 G772.G238 ADRP, PL2pro 

NSP4 G2362_Q2838 G2375.Q?85! 2390. ases 

NSP5 A2839. (3140 A2852_(3153 A2867 (03168 3CLpro 

NSP6 giu Qi) gis Qi» S3169 au 

NSP7 $3420. (3502 S353 355 S3448 (93530 

NSPS. Sa — Ag Bic que goa. Qs Primase 

NSP9 N3695. (93805 N3711_Q3818 N?726. 93833 

NSPIO A396 "0 A3319 (3953 3834.9 3968 

NSP11 oT? aL ae ae PO 3986 Short peptide at the end of ORFla 
NSPI2 — T3%!.Q4867 T?959. 880 173964. (34895 RdRp 

NSPI3. S#868.Q5464 A4881_Q5477 A896, 492 HEL, NTPase 

NSPl4  AS6.Q*»? A548. 03995 A593. georo ExoN, NMT 

NSPIS  S599.Q6*! S596 (96330 S6011_Q6348 NendoU 

NSP16  S92.y$95 S6331 K 6631 S6349 76646 OMT 


^ADRP ADP-ribose 1-phosphatase, PL2pro papain-like protease 2, 3CLpro coronavirus NSP5 protease, Hel helicase, 
NTPase nucleoside triphosphatase, ExoN exoribonuclease, NMT N7 methyltransferase, NendoU 
endoribonuclease, OMT 2’ O-methyltransferase 
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